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Recently, several groups have investigated quantum analogues of random walk algorithms, both 
on a line and on a circle. It has been found that the quantum versions have markedly different 
features to the classical versions. Namely, the variance on the line, and the mixing time on the 
circle increase quadratically faster in the quantum versions as compared to the classical versions. 
Here, we propose a scheme to implement the quantum random walk on a line and on a circle in an 
ion trap quantum computer. With current ion trap technology, the number of steps that could be 
experimentally implemented will be relatively small. However, we show how the enhanced features 
of these walks could be observed experimentally. In the limit of strong decoherence, the quantum 
random walk tends to the classical random walk. By measuring the degree to which the walk remains 
'quantum', this algorithm could serve as an important benchmarking protocol for ion trap quantum 
computers. 

PACS numbers: 03.67.Lx, 32.80.Pj 



I. INTRODUCTION 



The idea that a computational device based on the 
laws of quantum mechanics might be more powerful than 
a computational device based on classical mechanics has 
been around for about two decades 0]. The study of 
computational devices based upon quantum mechanics is 
known as quantum computation. For an introduction to 
the field, see for example Nielsen and Chuang ||^. Active 
research in this field has exploded since the discovery by 
Shor 1^ that a quantum computer could, in theory, factor 
large semi-primes exponentially faster than can currently 
be done on a classical computer. Since Shor's algorithm, 
Grover has devised an algorithm which can, in principle, 
search an unsorted database quadratically faster than 
any classical algorithm Q. However, new quantum al- 
gorithms which out perform their classical counterparts 
are proving difficult to find. One path which is being 
followed to find new quantum algorithms involves look- 
ing at effective classical algorithmic techniques, and try- 
ing to adapt them to quantum computation. Classically, 
the random walk has found applications in many fields 
including astronomy, solid state physics, polymer chem- 
istry and biology. For a review of the theory and appli- 
cations for random walks, see for example Barber and 
Ninham Q . The hope is that a quantum version of the 
random walk might lead to applications unavailable clas- 
sically. Quantum random walks have been investigated 
by a number of groups |^, ||, ||, ^ In this paper, 
we propose a scheme to implement the discrete quantum 
random walk on a line and on a circle using an 
ion trap quantum computer. For a review of ion trap 
quantum computation see Wineland et al. 1 12 



With 

current ion trap technologies, it will not be possible to 
implement a large number of steps in the walk, however 



it should be possible to implement enough steps to ex- 
perimentally highlight the differences between the classi- 
cal and quantum random walks, providing an important 
proof of principle. 

The structure of this paper is as follows. In Sec. II 
we review the simple models of random walks on both a 
line and a circle, highlighting the differences between the 



III 



classical and quantum versions in both cases. In Sec 
we discuss how we shall be representing the algorithms 
in an ion trap quantum computer. We then discuss the 
pulses required to evolve the system, first for the walk 
on the line, and then for the walk on a circle. Finally, 
m Sec. we discuss a relatively simple measurement 
procedure which can be used to highlight the difference 
between the classical and quantum random walks. 



II. CLASSICAL VERSUS QUANTUM RANDOM 
WALKS 

Classical random walks can take many different forms, 
starting from the simple discrete random walk on a line, 
to random walks on graphs, to continuous time random 
walks, such as brownian motion. In this paper we are 
only considering discrete time, discrete space, random 
walks on a line and on a circle. 



A. Classical walk on a line 

Imagine a person standing at the origin of a line with 
a coin in their hand. They fiip the coin, and if it comes 
up heads, they take a step to the right, if it is tails, they 
take a step to the left. They then repeat this procedure, 
hipping the coin, and taking a step based on the result. 
The probability, P/v {d) , of being in a position d after N 
steps is 
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PN{d) = ^ 



N 

d+N 



(1) 
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TABLE I: The probability of being found at position d after 
A'' steps of the classical random walk on the line. 



Table. || contains the probabilities for the first few val- 
ues of N. The non-zero elements of the distribution are 
simply terms from Pascal's triangle, divided by the ap- 
propriate factor of two. There are two features of this 
random walk that we would like to compare to the quan- 
tum analogue. Firstly, the mean of the walk is zero. This 
is intuitively obvious, we are using a fair coin, so we are 
as likely to step left as we are to step right. The other 
property of the distribution that we are interested in is 
the standard deviation. It is not hard to calculate that 
the standard deviation of this distribution, ac, is given 

by 



TV. 



B. Quantum walk on a line 



(2) 



Now let us consider a quantum version of the walk on a 
line. The first modification we can make is to replace the 
coin with a qubit. In this paper, we shall be representing 
the two levels of the qubit with the states | |) and | t) 
rather than |0) and |1). If we start with the qubit in the 
down state, and apply a Hadamard operation, we get an 
equal superposition of up and down, 
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If we were to measure the qubit, and step left or right 
depending upon the result, we would obtain exactly the 
classical walk described above. Now, rather than a per- 
son holding a coin, suppose we have a particle, whose 
motion is confined to one dimension. We can now treat 
the particle as a quantum system, and perform the quan- 
tum walk as follows. During each iteration, we apply the 
Hadamard operation, followed by the operation which 
steps right if the qubit is down, and steps left if the qubit 
is up. That is, we apply the operator. 



U 



(4) 



where p is the momentum operator of the particle con- 
fined to one dimension, and cr^ is the pauli-z operator 
acting on the qubit. Therefore, the state of the system 
after N steps is 
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TABLE II: The probability of being found at position d after 
N steps of the quantum random walk on the line, with the 
initial qubit state ^| j) -I- -^| j). 



where |^'o) is the initial state of the system. The mean of 
the distribution produced by this quantum random walk 
is not necessarily zero. It is dependent upon the initial 
state of the qubit. For example, choosing the initial state 
of the qubit to be down gives a non-zero mean after the 
second step. For the remainder of this paper, we shall 
only be considering the distribution created with the ini- 
tial qubit state -^| i) + "^1 T) which has a mean of zero 
for all values of N , 



l*o) = -^|0)(U)+*|t)). 



(6) 



Table. || contains the probability distribution associated 
with the first few states |^jv)- The non-zero elements of 
the distribution are no longer simply terms from Pascal's 
triangle which arose in the classical case. The deviations 
from the classical distribution are caused by quantum 
interference effects. Now it is no longer simple to calcu- 
late the standard deviation of the distribution. However, 
numerical simulations reveal that the standard deviation, 
CTq, is almost independent of the initial state of the qubit, 
and is approximately linear in N , 



-N. 



(7) 



The standard deviation is plotted in Fig. |l| up to = 40 
for both the classical and quantum walk distributions. 
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FIG. 1: Standard deviation for both the quantum and classi- 
cal random walks up to A = 40. 
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Clearly, the standard deviation is significantly different 
for the quantum and classical random walks on a line. 
Now let's consider the random walks which arise when 
periodic boundary conditions are applied to the random 
walks. 



C. Classical walk on a circle 

In the paper by Aharonov et al. they consider 
random walks on the circle, where the step size is an ir- 
rational multiple of tt. Here, we shall only be considering 
the simple distribution which arises when the step size 
is taken to be 7r/2. Let us assume that the particle is 
initially found, with probability one, at some point on a 
circle denoted by 6* = 0, 



Po{e = o) = 1. 



(8) 



After one step of the algorithm, the classical distribution 
is given by 



and after the second step, 

P2{0) = I 



= 0,7r 

2 



= 0,7r 

^ 2 



(9) 



(10) 



It is not difficult to see that the probability distribution 
for all subsequent odd number of steps will be given by 
Eq. (^) , and the distribution for all subsequent even num- 
ber of steps will be given by Eq. ( p^ ) . 



D. Quantum walk on a circle 

Let us consider the quantum random walk on a circle. 
Once again, we start with the particle at some point on 
a circle denoted hy 9 — 0, thus the initial probability 
distribution is given by Eq. (^. The probability distri- 
butions after one and two steps are also given by Eqs. (||) 
and ( p^ respectively, however after the third step, inter- 
ference effects results in the distribution 



Ps(0 = ^) 



1. 



(11) 



Calculation of the states after subsequent steps reveals 
that the quantum random walk around the circle, with 
a step size of 7r/2 is periodic with a period of eight. The 
eight probability distributions which arise are given in 
Table, ffli. 



III. 



IMPLEMENTING THE WALKS IN AN ION 
TRAP 



The analysis thus far has assumed that all operations 
can be applied without error and the particle can ex- 
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TABLE III: The probability of being found at position i 
steps of the quantum random walk on the cirle. 
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ist in position eigenstates. Now we shall relax these as- 
sumptions, and describe how the algorithm can be im- 
plemented in an ion trap. 

The ion trap provides a convenient setting for the 
quantum random walks we have described, as it contains 
the required discrete and continuous quantum variables. 
For the remainder of this paper we shall be discussing 
implementations based on a single ^Be"*" ion, confined in 
a coaxial-resonator radio frequency (RF)-ion trap, as de- 
scribed in and references therein. 

The preparation involves laser-cooling the ion to the 
motional and electronic ground state, |0)| |), as described 
in |l^ . A sequence of four Raman beam pulses are then 
applied to create the state (|a)| i) + | — Q;)| t))/\/2, 
where la) denotes the coherent state of the the oscillator. 



da;e^*"^="e-^(^-^"«)'|x) (12) 



ri/4 



and a = an + iai. 

The first pulse is a 7r/2-pulse which creates an equal 
superposition of |0)| |) and |0)| t)- A displacement beam 
is then applied which excites the motion correlated to the 
I I) internal state. The third pulse is a 7r-pulse which ex- 
changes the internal states, and finally the displacement 
beam is applied again. The combined action of the four 
pulses is to effectively perform the operator U, defined 
in Eq. (^. The quantum random walk on the line is ac- 
complished by repeating this sequence of pulses N times. 
Fig. 1^ contains the Wigner function obtained by tracing 
over the internal degree of freedom after five steps of the 
quantum random walk algorithm. 

The quantum random walk on the circle can be imple- 
mented in an ion trap by 'walking' the particle around a 
circle in phase space, rather than a circle in real space. 
In order to accomplish this task, we need to generate an 
operator of the form 



W 



(13) 



where a and correspond to the annihilation and cre- 
ation operators of the harmonic oscillator. This operator 
can be produced in an ion trap by applying far-detuned 
laser pulses to the ion |l5|, followed by a 7r/2-pulse. 
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Quantum Random Walk on a Line Probability of measuring ground state (Step = 2) 




FIG. 2: Wigner function of the particle after five steps of the 

quantum random walk on the line. (The electronic level of ^ 

the ion has been traced over.) 10 20 30 

time (units of t / co) 



IV. MEASURING THE WALKS 

Using current ion trap technologies, wave packet dis- 
persion is negligible , so the main source of decoher- 
ence is related to the internal levels of the ion. Decoher- 
ence of the electronic levels of the ion during the appli- 
cation of the algorithm has the effect of gradually trans- 
forming the quantum random walk to the classical ran- 
dom walk. Rather than considering this to be a negative 
effect, we can measure the degree to which the ion is act- 
ing as a quantum variable rather than a classical variable, 
and thereby effectively measure the level of decoherence 
in the ion trap. 

The scheme that we envisage for measuring the ran- 
dom walk utilises similar operators to those employed in 
the application of the algorithm. After applying the ran- 
dom walk sequence for some number of steps, the internal 
state of the ion is decoupled from the motional state by 
an appropriate Raman pulse. An effective operator such 
as exp{ipay) is applied, before finally measuring the in- 
ternal state of the ion. Thus we are using the internal 
state of the ion to supply as with information about the 
motional state. 

In the case of the walk on the line, suppose we decouple 
the internal state from the motional state by measuring 
whether the ion is in the state | t) or | 4). We then apply 
the operator 

M± = e^ip<^y. (14) 

The positive Hamiltonian is applied upon obtaining the 
results I t)j whilst the negative Hamiltonian is applied 
otherwise. Finally, we again measure the internal state of 
the ion. If the quantum random walk has experienced no 
decoherence, then we measure | |) with the probabilities 
given by the solid line in Fig. ^, whereas if the ion suffers 
complete decoherence we would expect to measure | J,) 
with probability of one half. 

A similar scheme can be used to measure the level of 
decoherence in the quantum random walk on the circle. 
Fig. H again depicts the probability of measuring the ion 



FIG. 3: Probability of measuring the ion in the ground state 
after applying the random walk for a time t/u, decoupling the 
internal and motional states, and applying the measurement 
operator Af . 

in the ground state after decoupling the internal and mo- 
tional states, however this time we then apply the oper- 
ator 

D = e"*«. (15) 

In this case, because we have total destructive interfer- 
ence of certain paths during the walk, the deviation of 
the quantum to classical walk is much larger at certain 
stages of the walk. 
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FIG. 4: Probability of measuring the ion in the ground state 
after applying the random walk on a circle for a time irt/uj, 
decoupling the internal and motional states, and applying the 
measurement operator D. 
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V. DISCUSSION 

We have described ion trap implementation schemes 
for quantum random walks, both on the lino and on the 
circle. We have also suggested a measurement process 
which allows the enhanced features of these walks to be 
experimentally observed. 

At this point, it is unclear whether quantum random 
walks will have any useful algorithmic applications. How- 
ever, we believe that they can provide a benchmarking 



protocol for ion trap quantum computers, and perhaps 
other implementation schemes which combine continuous 
and discrete quantum variables. 
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